Prepared for submission to JHEP 



Conductivity and quasinormal modes in holographic 
theories 



M.A. Stephanov, Y. Yin 

Department of Physics, University of Illinois, Chicago, IL 60607-7059, USA 
E-mail: misha@uic.edu, yyin3@uic.edu 

Abstract: We show that in field theories with a holographic dual the retarded Green's 
function of a conserved current can be represented as a convergent sum over the quasinormal 
modes. We find that the zero-frequency conductivity is related to the sum over quasinormal 
modes and their high-frequency asymptotics via a sum rule. We derive the asymptotics 
of the quasinormal mode frequencies and their residues using the phase-integral (WKB) 
approach and provide analytic insight into the existing numerical observations concerning 
the asymptotic behavior of the spectral densities. 
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1 Introduction and summary of results 

The transport properties of the strongly coupled quark-gluon plasma (sQGP) created at 
RHIC [1-4] attracted much attention recently. One of the most important transport pa- 
rameters is the conductivity a associated with a conserved vector current. For example, the 
quark current conductivity is an indicator of deconfinement. Furthermore, the conductivity 
of the current of light quarks can be related, via Kubo formula, to the soft limit of the 
thermal photon production rate by QCD plasma. In addition, the Einstein relation equates 
conductivity to the product of quark susceptibility and quark diffusion constant. For heavy 
quarks, the diffusion constant is an important quantity characterizing medium effects on 
quark propagation. The results of application of various phenomenological models to that 
problem is best expressed in terms of the diffusion constant [5]. 
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Due to strong coupling, the calculation of conductivity in QCD at temperatures rele- 
vant to the experiments is a challenging task. Lattice calculations, being restricted to the 
finite interval of Euclidean time, require analytic continuation to infinitely large real time 
in order to determine transport coefficients such as conductivity. Some interesting results 
have been obtained assuming certain analytic behavior [6]. Clearly, it would be greatly 
helpful to better understand the analytic properties of the current-current correlator and 
find generic, model-independent constraints such as sum rules on the conductivity. 

As a step towards this goal, in this paper, we consider the class of quantum field 
theories in 3 + 1 space-time dimensions whose correlation functions can be computed via 
AdS/CFT duality [7-9]. Such theories have been employed to describe thermodynamics 
and transport in strongly coupled regime of QCD. We consider a generic gravity dual set- 
up satisfying some mild technical assumptions as considered in Ref. [10]. One can show 
that the retarded Green's function at vanishing momentum Gji{u}) of a conserved vector 
current calculated from such gravity background is a meromorphic function with infinite 
number of simple poles located in the lower half-plane [10]. In the context of gauge/gravity 
correspondence, those poles are referred to as quasi-normal modes [11-16]. 

We show that the Green's function Gr{uj) can be represented as a convergent sum over 
its poles: 

+ 1 . (1.1) 

in terms of conductivity a and the residues r„. The real coefficient C depends on the 
definition of Gf({uj), however, its temperature-dependent part is fixed by parameters of the 
second order hydrodyanmics [17] (and equals a times tj defined in Ref. [18]). Since Gr{oj) 
has a "mirror" symmetry: G^(ct;) = Gr{Gj) where Co = — w*, if a;„ is a pole of G/j(ct;), so 
is Cjn-^ The index n in the infinite sum in Eq. (1.1) counts poles located in the fourth 
quadrant of the complex oj plane. The contribution of the mirror poles is the second term 
in the infinite sum. 

Expansion of Gr{uj) in terms of poles corresponding to quasinormal modes has been 
suggested in Ref. [19, 20]. However, to avoid ambiguities, one needs to show that the 
summation in the expansion is convergent for any finite u away from the poles. To this 
end, we have established the convergence by determining the large n asymptotics of both 
LOn and r„ (by extending previous work [21] on asymptotics of a;„): 

LOn — ^ nLOQ + A., r„ — )• Kloq, when n — )• oo. (1.2) 

The complex numbers loq and A are sometimes called (asymptotic) "gap" and "offset" 
of the quasi- normal modes respectively [21]. The coefficient K is related to the leading 
asymptotic behavior of Gr, which in the deep Euclidean regime co — )■ too is given by the 
operator product expansion (OPE): 

GR{iuJE) ^ -2KujElogUJE, UJE ^oo. (1.3) 

^For simplicity, we assume there is no pole located at the negative imaginary axis thus w„ Cun- If there 
is, the modification of our treatment here is trivial. 



Gr{uj) = -iauj + Cuj'^ + ^^yZ 
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The constant K is proportional to the number of the charge carrying degrees of freedom. 

Finally, by matching the asymptotic behavior in the deep Euclidean regime of the 
representation (1.1) to the OPE (1.3) we derive a relationship between the conductivity 
and the quasinormal modes: 

a = -Kim (wo + 2A) + 2 ^ Im (r„ - Kujq) . (1.4) 

n 

This paper is organized as follows. Sec. 2 presents the derivation of the representation and 
the sum rule. In Sec. 3, we establish the asymptotics of a;„, r„ using the WKB approxima- 
tion. In Sec. 4, we investigate how the sum rule is saturated by studying the "soft-wall" 
model [22] at finite temperature numerically. We summarize and explain qualitatively and 
quantitatively how the asymptotic behavior of Un, Vn is related to the "damped oscillating" 
behavior [19] of spectral densities in Sec. 5. In Appendix A, we clarify a subtle point in the 
holographic calculation of the retarded correlators in the lower half of the complex u plane. 
In Appendix B we derive the Stokes constant formula we used in the WKB calculation. 
We also formulate a family of f-sum rules from holography in Appendix C. 



2 The derivation of the representation and the sum rule. 
2.1 The representation 

We study the retarded Green's function Gfj,{uj) of a spatial conserved vector current oper- 
ator at zero three-momentum and the corresponding spectral function p{u): 

Giiiuj) = -i [ dte'^H{t){[j\t),j\{))\), P{^) = -lmG^j(w). (2.1) 



We assume that the quantum field theory under consideration has a holographically dual 
description. As discussed in Ref. [10], Gr{uj) calculated from holography, is a meromorphic 
function on general grounds. We could thus consider the following Mittag-Leffler expansion 
of Gfi{uj) modulo contact terms: 



Gr = — 2— = -i-+ujy 



+ — 



U)n (w - UJn) W„ {UJ - U)n) 



+ P{u), (2.2) 



where -P(w) is a polynomial of u. The scaled retarded Green's function Gr{uj) defined in 
Eq. (2.2) has a pole at w = with residue related to the conductivity by the usual Kubo 
formula: 

a = hm ^ (2.3) 

ti;->-0 LO 

while the quasinormal mode residues are defined as 

r„ = lim {uj - ujn)GR{uj). (2.4) 

In the deep Euclidean region, in accordance with the operator product expansion (OPE), 
Gr{oj) has the following asymptotics: 

G'ij(ia;£;) = 2i^ log a;^; -|- const -|- 0(ti;^^), w^; — )■ oo (2.5) 
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where the leading contribution is from the unit operator. Here we have used the relation 
GR{ioJE) = —Ge{oje > 0) where Ge{oje) is the Euclidean correlator ^ and the OPE of 
GEioJE)- When writing down Eq. (2.5), we have assumed that the lowest dimension of 
those non-trivial operators entering the OPE of Ge{oje) is no less than 2. That fact is 
crucial for subsequent discussions. 

Since P{oj) is a polynomial, the logarithmic behavior in Eq. (2.5) should be matched 
by the summation over pole contributions in the representation (2.2). Thus the number 
of poles has to be infinite. As we will show in the next section, uon, fn have the asymptotic 
behavior given by Eq. (1.2). As a result, for any finite lo away from uJn{ojn), the summation 
of r.„/a;„(a; — ujn) in Eq. (2.2) is convergent. 

We also note from Eq. (2.5) that the polynomial P{uj) cannot grow faster than a 
constant. Consequently, it should be a real constant, i.e. C, as required by the "mirror" 
symmetry of Gr{uj). We thus establish the representation (1.1). 

2.2 The sum rule 

In order to derive the sum rule relating conductivity a to the quasinormal modes, we shall 
match the asymptotic behavior of representation in Eq. (2.2) to the OPE Eq. (2.5). 
To facilitate the matching, we apply "Borel" transformation [23] defined by: 

l^n / d \^ UlE 

^iAb = (^7^ j > when o;^ ^ +00, n ^ oo, ^ = l/t^, (2.6) 

to Gfj,{uj) in the deep Euclidean region: 

Bi/t^GR{iujE) = -tBCT-itBY.{rne-"^-'^ +fne-'''"'^) = -tBa+2tBj2^m (r^e"^'^"*^) . 

n n 

(2.7) 

All relevant formulas for the Borel transformation are listed in Appendix. C. For any 
positive tB, the sum in Eq. (2.7) is convergent since Ima;„ < 0. Applying the Borel 
transformation to the asymptotic expansion (2.5), we obtain for small tB ■ 

Byt^GR{iL0E) = -2K + O{tl), when ^ 0+. (2.8) 

Matching Eq. (2.7) and Eq. (2.8) at small tB, we find: 

K 



(7 = 2 lim 



Im rne~*'^"*s + 



tB 



(2.9) 



One can check, using Eq. (1.2), that when tB — )■ 0"*", the l/t_B divergence in Eq. (2.9) is 
canceled as it should be. 

Using the asymptotic behavior of r„ in Eq. (1.2) we can evaluate R.H.S of Eq. (2.9) 
by rearranging the infinite sum as 



£7 = 2 lim 



Im y (r„ - Kojo) e'"'^"^^ + A'lm V woe-^'^"*^ + — 



(2.10) 



^We analytically continue the Euclidean correlator Gsii^E) from the discrete set of Matsubara frequen- 
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The summation of (r„ — Kojq) is convergent due to Eq. (1.2) (see discussion in Sec. (3)). 
Consequently one can exchange the sequence of summation and taking — ?• 0"*" hmit. The 
second sum in Eq. (2.10) can be evaluated explicitly for Ib — >• 0"*". Its divergence —K/ts is 
cancelled by the last term in Eq. (2.9) and the remaining finite part can be obtained using 
asymptotics of uj^'- 



lim 

ts->0+ 



Im y (Joe-''^"*^ + — 

^ tB 



lim 



oo ^ 

Im y L^oe"'("'^°+^^*^ + — 

^ tB 



n=l 



lim 



Im 



Woe 



-iAt 



B 1 



(2.11) 



Expanding the expression in brakets around Ib = 0^, taking the limit and substituting into 
Eq. (2.10), we obtain the sum rule Eq. (1.4). 

2.3 J\f = 4 SYM theory in large Nc, strongly coupling limit as an example 

In Ref. [24], Gji{u}) in A/" = 4 SYM theory is derived in large Nc, strong coupling limit 
using AdS/CFT correspondence: 



\ 2-kT ^ {2-kTY 



UJ 



(1 



A-kT 



+ '4> 



AttT 



(2.12) 



with ip the logarithmic derivative of the gamma function. This retarded correlator has a 
quasi-normal spectrum with: 



UJr. 



2ttT{1 



in. 



iV2 

167r2 



(1 - ^)T. 



(2.13) 



Therefore, for this theory, ujq = 2(1 - i)TTT, A = and K = N^/32tt^. From Eq. (2.12), we 
also have a = A'^T/lGvr. One sees immediately that sum rule (1.4) holds. 



3 The asymptotics of quasinormal frequencies and residues 

3.1 The Green function in the holographically dual description 

To calculate Gr{uj) using gauge- gravity/holographic correspondence we need to consider 
the second order variation of the 5-dimensional bulk action with respect to the bulk gauge 
field Vm dual to the vector current in the boundary theory. The relevant part of the 
bulk action has the usual Maxwell form: 

S = J d^x^e-'t'VMNV^^ (3.1) 

where g'^ is the 5D gauge coupling, cf) is the background scalar field which, in general, is a 
combination of dilaton and/or tachyon fields, corresponding to the conformal and/or chiral 

^As a side, we note that the radius of convergence of the Taylor expansion in ts is |27r/(^o| because of 
a pole at = 2tt/luo- That suggests that |a;o|/27r sets a scale below which the asymptotic expansion (or 
OPE) of GR{iOE) will be broken. 
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symmetry breaking, and Vmn = QmVn — OnVm- We consider the most general metric 
(up to general coordinate transformations) possessing three-dimensional (3D) Euclidean 
isometry: 

= e2^(^) {hdf - dx^ - h-'dz^) . (3.2) 
The equation of motion resulting from the action (3.1) reads: 

d,{he^d,V) + uj^h'^e^V = (3.3) 

where V = V\ and B = A — ^. As usual, the thermal bath is represented by the black 
brane, corresponding to a real positive zero of /i(z), with temperature T given by:^ 

Here and hereafter a prime denotes the derivative with respect to z. Also in this section, 
we will set ttT = 1 for convenience, i.e., we will measure all dimensional quantities in units 
of ttT. We require the background to be AdS in the asymptotic at the boundary: 

i?(z) = — logz, when z — )• 0. (3.5) 

Then z = and z = zh are two regular singular points of Eq. (3.3). The retarded Green's 
function, up to a contact term, is given by the standard holographic prescription [25, 26]: 



Gr{uj) = lim 
55 -'^0 



he — + cj log z 



(3.6) 



where V- denotes the Frobenius power series solution near zh of indicial exponents —ioj/A'kT, 
i.e., ~ (z — z//)~*^/^'^-^[l + 0{z — Zh)]-, corresponding to an in-falling wave. Further 
assuming the Frobenius power series solutions at z = and z = zh have an overlapping 
region of validity along the real axis Q < z < zh, one can show, along the lines of Ref. [10], 
that Gr{uj) is a meromorphic function^. 

3.2 Near-boundary asymptotics 

To establish the large n asymptotics of quasinormal mode parameters in Eq. (1.2) we first 
introduce the book-keeping Schrodinger coordinate 

and a wave-function-like 

^r(^) = e^(")/V(z). (3.8) 
Then Eq. (3.3) is brought into the standard Schrodinger-like form: 

d^^ii) 



de 



+ {u' - U{0) ^(0 = (3.9) 



■'For simplicity, we assume that h{z) has only one real positive zero. 

^That conclusion from Ref. [10] has some uncertainties at a set of discrete frequencies uj = —2in-KT 
where the difference between two indicial exponents r+ — r_ is an integer. We will settle that subtle issue 
in Appendix. A 
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where the potential U{z), as a function of z, is given by 

U{z) = 



B'V B" h'B' 
+ — + 



2 J 2 2h 
Near the boundary z = ^ = 0, Eq. (3.9) becomes 



(3.10) 




■ y-- - ^ j m = 0' (3-11) 
where 1^0 = 1 (the same as the spin of the fluctuations we are studying).^ Its solutions are: 



m = A^{u)f-fH^^){u:0 + A^{u:)f-^H'^}{uO (3.12) 

with H^^ , hIj''^ denoting the Hankel functions of the first and second kind respectively. 
Using the definition (3.6), we can calculate the asymptotic behavior of Gr{ijj) from the 
solution (3.12): 



Gr{u:) = 2K 



7 + log(a;/2) + 



+ 0{uj-'^), with K = ^. (3.13) 



2tan[a;P(tj)] 

Here, 7 is the Euler-Mascheroni constant. The function T>[uj) is defined by 



1 

27i 



-2iujV{uj) 



(3.14) 



It will be determined by applying the in-falling wave boundary condition near the horizon. 
The poles of Gr{uj), ujn, as well as the residues r^, for sufficiently large n can be determined 
from Eq. (3.13): 

ttK 

UJnV{u:n) = rnr, rn = ^ , , . (3.15) 

One may note that corrections in Eq. (3.13) are 0(a;~^) as we are required to match 
Eq. (3.13) with the OPE results Eq. (2.5)^. As a result, a;„,r„ calculated from Eq. (3.15) 
are accurate up to (including) the order of relative to the corresponding leading large 
n results. 

To determine uin and r„ from Eq. (3.15) for large n, we need to know asymptotics of 
P(a;). We shall determine it by following the solution of the Schrodinger equation along 
a path from z = to z = zh where we apply the in-falling wave boundary condition. We 
shall use the WKB solution along that path. Thus it is important that the region, which we 
denote Rq, where \z\ <^ 1, does overlap with the region Ri, defined by \io'^\ |f/(Oli where 
WKB approximation is applicable, for sufficiently large oj. We denote that overlapping 
region by i?2- In R2 ^ 1. (The above definitions of the regions are summarized in 
Table. 1). 



"^Our formalism below can be readily generalized to other values of vq- 

^From the gravity side, that condition will be satisfied if L''(^) — (fQ — 1/4)/^^ is bounded near the 
boundary. 
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Table 1. The definition of different regions. 



Rq 


z <C 1 and Eq. (3.9) is well approximated by Eq. (3.16). 


Ri 


1^^ ^ f^(C) ^^(^ one can use the WKB approximation. 


R2 


Rq n Ri where we will match Eq. (3.21) and Eq. (3.17). 



Finally, due to the asymptotic behavior of the Hankel functions at large argument, 
hI^\x) K \ —exp(ix -i6), hI^^ \ — exp(-ix + i6) (|x| > 1), (3.16) 

V VrX \ TTX 

where 5 = {2u + l)'7r/4, we find for large |a;^| (i.e., in the region i?2): 

^(0 = A+e^'^«-*'^o + A_e-^^«+^''o, with e G i?2 (3.17) 
where 5q = Stt/A. We shall match the WKB solution to this asymptotics. 
3.3 The WKB approximation 

In region Ri one can use WKB approximation (also known as the phase integral method [27]) 
to solve Eq. (3.9). The application of the method to calculating the asymptotic quasi- 
normal modes is reviewed in Ref. [21]. 

The two linearly independent WKB solutions are given by Q~^^'^ exp[ibz JdS^Q], where 

Q\^) = u:^-U{0 (3.18) 

These solutions are singular at points where Q = - the turning points. At a turning 
point the WKB approximation breaks down. We shall assume that a generic Schrodinger 
potential grows as U{z) ~ when |z| is large, where m is a positive real number. Then in 
the limit of large |a;| , there will be turning points determined by the condition — z™ = 0, 
which has multiple solutions. We denote one such turning point by zt and map it to the 
Schrodinger coordinate S.T = Ci^r)- 

The exponential in the WKB solutions is purely oscillatory along an integral curve 
defined by condition 

Im [ d^'QiC') = (3.19) 

This condition defines the anti-Stokes line(s), AS, with respect to the point (,TizT)- For a 
simple zero ^t{zt) of Q{^{z)), there will be three anti-Stokes lines ASi,AS2, AS3 emanating 
from it. We choose zt (or ^t) among solutions of Q{S,{z)) = by requiring that ASi has 
an overlap with R2 while AS2 ends on zh as illustrated by Fig. 1(a). As it will become 
clear soon, for large the existence of such a turning point zt (or ^t) is necessary to 
solve the condition (3.15). 

We shall define point ^00 as the limit 

eoo= lim eT = ,hm / -— (3.20) 
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(a) (b) 

Figure 1. The Stokes diagram for uj = 10(1 — i)TrT of the pure AdS black-hole background in the 
complex z plane (left panel) and in the complex ^ plane (right panel). We set ttT = 1. Anti-Stokes 
lines ASi,AS2,AS3 are blue solid lines while Stokes lines 31,82,83 are red dashed lines. The path 
r connecting ASi and AS2 is plotted as a thick black line. To determine the asymptotics (1.2), 
we connect the region R2 (around 2 = 0) and horizon zh via the path along ASi, T and AS2- We 
sketch the boundary of Ri (where the WKB approximation applies) schematically using the green 
thick dashed lines. For completeness, we plot Stokes lines and anti-Stokes lines emanating from 
another turning point z!j, next to zt- 



which exists if the integration in the above equation is convergent. To avoid ambiguity, we 
specify the path of the integration in Eq. (3.20) to be ^5*1. Since, in the large |a;| limit, the 
region Ri extends to the origin, the origin z = and the turning point zt are connected 
via ASi in that limit. 

We now start using the WKB approximation to solve Eq. (3.9). Along ^l^i, we can 
express ^'(^) using the standard WKB approximation: 

q,^^^) = 4^e*/L '^'Q + ^e-'^L -^^'^ « J^e^^^"?-) + 4^e"-«-«-) (3.21) 

vQ VQ VQ VQ 



as long as ^ stays in region Ri. For convenience, we choose .^00 to be the lower limit of 
the "phase integration" in Eq. (3.21). A different choice of the lower limit of the integra- 
tion would not affect our final results, but would complicate their derivation. Matching 
Eq. (3.21) with Eq. (3.17) in region R2, we have: 

.A B 

Similarly, along another anti-Stokes line AS2 in region Ri, we could write ^(^) as: 

^^(^) ~ ^e^'^(«-«°°) + -^e-^'^(«-«°°). (3.23) 

vQ VQ 



As both ^siO and ^ciO represent the same solution of the Schrodinger equation (3.9) in 
different Stokes domains, C± can be expressed as a linear combination of B±. To see that, 
we trace ^siO along a path F connecting ASi and AS2 while staying in Ri as illustrated 
in Fig. 1. Along F, the "in-falling wave" term e^^^^~^°°^ is (exponentially) dominant over 
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the second term e ^'^^^ Therefore 5+ must not change along F, no matter the value 
of 5„: 

C+ = B+. (3.24) 

On the other hand, if -B+ = 0, then i?_ cannot change along F, i.e., C_ = if i?+ = 0. 
As a result, we can express C_ as: 

C- = B^ + SB+ (3.25) 

where the multiplier S is the Stokes constant [27] with respect to point .^oo- The phe- 
nomenon that the coefficient of the subdominant solution is shifted by a product of S and 
the coefficient of the (unchanged) dominant term is the well-known Stokes phenomenon.® 
In addition, near zh, we have from Eq. (3.7): ^ ~ — log{zH — z)/4:. Then selecting the 
solution V-{z) in Eq. (3.6) is equivalent to imposing the infalling wave condition: 

C_ = (3.26) 

along AS2 [21]. From Eq. (3.25), we obtain: 

^ = -4 = e"'°^^"'''. (3.27) 
B- S 

Substituting the above equation (3.27) in Eq. (3.22) and using the definition (3.14), we 
establish an asymptotic expression for P(a;): 

ojV{oj) = - ^ (2z/o - 1) - ^ log S (3.28) 



and from Eq. (3.15) the asymptotic behavior of ujn and r^: 

Wo, rn = K iuJo + 



n + -^log5 + ^(2z/o- 1) 
zvr 4 



icon dS 



2TTSduj 



(3.29) 



where 



wo = vr/^oo- (3.30) 



3.4 The Stokes constant 



To gain insight into how and whether S should (or should not) depend on w, it is useful to 
think of the Stokes phenomenon in the following way [28]. Both WKB solutions Eqs. (3.21) 
and (3.23), describe the same exact solution of the Schrodinger equation in two different 
Stokes sectors around the turning point. These approximate solutions are multivalued 
functions with a branching singularity at the turning point, Q{S,t) = 0. However, the 
exact solution of the Schrodinger equation is analytic at a regular point, such as the turning 
point. In order to match the absense of the branching singularity in the exact solution, 
the WKB solutions must compensate their discontinuity along a path winding around the 
turning point by a corresponding discontinuity in the coefficients B and C. This is the 



*This shift occurs along F discontinuously at the crossing of tlie Stokes hne separating the Stokes 
domains. 
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essense of the Stokes phenomenon. For an isolated regular turning point this argument 
leads to the well-known value oi S = i. 

More importantly, this argument sheds light on the reason why the Stokes constant 
should have a different value in a special case when the turning point approaches a singu- 
larity of the Schrodinger potential in the limit \uj\ — )• oo, as it does in our case. Since, in 
this case, winding around the turning point, while staying in Ri, requires winding around 
the singularity also (as well as other turning points). If the singularity is a branching point 
of the exact solution, the discontinuity across the cut is reflected in the value of the Stokes 
constant, which thus depends on the nature of the singularity. 

In fact, since we assume that, at large z, [/ ~ and z/h vanishes to guarantee 
the convergence of the integration in Eq. (3.20), ^oo will always be a singular point of 
Eq. (3.9). If ^oo is a regular singular point of Eq. (3.9), ^'(^) can be expressed as a linear 
combination of two Frobenius series solutions: (^ — ^oo)'^'''(l + — Coo)) with f± being 
the indicial exponents and f+ + f- = 1. Taking the WKB solution around the point 
^oo and matching the discontinuity of the exact solution, one finds the Stokes constant: 
S = 2i cos[7r(/_|. — /_)/2] [28] as we explain in detail in Appendix B. If the indicial 
exponents f± are independent of w, the resulting Stokes constant has no uj dependence 
either. 

Even if ^oo is an irregular singular point of Eq. (3.9) one may still expect that deter- 
mining S, though more involved, is still possible, perhaps along the lines of Ref. [29, 30] 
(see also Appendix B). From a more practical point of view, which we take in Sec. 4, even 
if one has not found an easy way to determine S in that situation, one could attempt to fit 
asymptotic behavior of Un numerically using Eq. (3.32). If the quality of the fit is good and 
the resulting loq is close to the analytical expectation given by Eq. (3.20), then it is very 
likely that S will approach a constant in large \uj\ limit. In fact, that is what we observe 
for the soft- wall model at finite temperature (see Sec. (4) below). 

In conclusion, we anticipate, on general grounds, that for a large class of theories the 
Stokes constant S" is a finite constant in the large |a;| limit. As a result, the asymptotics (1.2) 
are established. 

Furthermore, to show that the summation over (r„ — Kuq) is convergent in the sum 
rule (1.4), we need to show that summation over S^^dS/du) terms is convergent. For 
sufficiently large n, we can replace the summation with integration. Then the existence 
of the large |w| limit of log 5" would imply that the integration of S~^dS /doj is convergent 
and complete the derivation of the conductivity sum rule (1.4). 

3.5 Examples and comparisons 

The authors of Ref. [21] have considered the cases that the Schrodinger Equation (3.9) 
takes the form: 

^r(^) = 0, when \z\ oo. (3.31) 



+ 



UJ 



(^-^oc)2 
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Then f± = ±1^00 — 1/2 and we have S = 2z cos(7ri/oo)(see also Ref. [27]). Consequently, we 



read from Eq. (3.29) that: 



A 



i 1 

— log(2zcos(7ri/oo)) + -(2z>o - 1) 
zvr 4 



(3.32) 



in complete agreement of the results of Ref. [21] obtained by using the properties of the 
Bessel functions^. For that reason, the first part of Eq. (3.29) is a generalization of the 
previous work. Although asymptotic behavior of r„ can be calculated straightforwardly 
from the WKB approximations, the expression of r„ in second part of Eq. (3.29), to the 
best of our knowledge, is new. 

Finally, let us check the results obtained in this section in the case of pure thermal AdS 
black-hole background. In that case, B{z) = — log z, U (z) « — 5z^/4 and ^ « + l/(3z^) 
when \z\ is large. In that limit, the Schrodinger equation (3.9) is reduced to Eq. (3.31) 
with = —1/3. From Eq. (3.32) and recalling that vq = 1, we obtain A = due to the 
cancellation between two terms in Eq. (3.32). Moreover, if h{z) is a polynomial, as it is in 
the case at hand, h{z) = 1 — z^, there is a simple way to evaluate ^00 in Eq. (3.20): 



Co 



Zoo 

-^oo Hz) 



1 

2jc 



dz 



1 



h{z) 



-TTl 



E 



1 



h'izh] 



(3.33) 



In the first equality we have used the property h{z) = h{—z). The contour C is chosen 
to connect ibzoo by a straight line and a large semi-circle centered at the origin of the 
complex z plane. As [zool 00 when |a;| — t- 00, the contribution from the integration 
along the semi-circle vanishes for h{z) being a polynomial of z^. Applying the Cauchy 
integral theorem to the integral, we obtain the rightmost expression in Eq. (3.33). The 
summation here denotes the summation over all z^s, the zeros of h{z), enclosed in the 
contour C. In particular, for the pure AdS black-hole background, ArgZoo = 7r/12 as can 
be seen in Fig. 1(a). Therefore z = are the zeros of h{z) enclosed in the contour C. 

Consequently, ^00 = (1 + ^)/2^ thus ujq = vr/^oo = 2(1 — i)-KT (we restored the units which 
were set by irT = 1 in this Section). One can check that ujn,rn given by Eq. (1.2) coincide 
with the quasi-normal spectrum of A/" = 4 SYM theory given by Eq. (2.13). 



4 Examining the sum rule in the "soft-wall" model at finite temperature. 

In this section, we will examine the sum rule (1.4) with the "soft-wall" model [22], a 
holographic QCD model, at finite temperature. With (/)(z) = cz^ and A{z) = — log z, the 
"soft-wall" model reproduces the Regge-like trajectory of the vector mesons, = 4nc, at 
zero temperature [22]. Studying that model at finite temperature can provide a non-trivial 
check of the sum rule (1.4). It can also illustrate how the dissociation or "melting" of the 
bound-states is related to the increase in conductivity, the phenomenon which is relevant 
to the transport properties of sQGP. 

Dimensionless ratios of physical quantities in the "soft-wall" model at finite temper- 
ature are fully controlled by the dimensionless parameter c = c/{'irT)'^. To make the 

®We have converted the results of Ref. [21] into the notations used in this paper. 
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0.25 0.30 0.35 0.40 0.45 0.50 

T(GeV) 



Figure 2. a{Green), s(dashed), si(blue), and S2(red) divided by the constant 2TrK as a function 
of T. To evaluate S2, we truncated the summation at n = 5. We also plot the contribution from 
the A term to the total sum normalized by 2TrK: sa = —2KlmA in black. 

connection with the real world more tangible, we set the overall scale by taking c = 2.54 
GeV to fit the mass of the J/ip at zero temperature. Such a choice has been used in 
Ref. [31, 32] to study the thermal charmonium spectral functions. 

Following Ref. [31-34]), we assume pure black-hole metric background, i.e, h{z) = 
1 — (zttT)^. We have calculated the first five quasi- normal modes ujn and the corresponding 
rescaled residues r„ numerically for temperature T between 250 MeV and 500 MeV. For 
the soft-wall model at finite temperature, the point is an irregular singular point of 
Eq. (3.9). As we have explained in the previous section, we fit ujn using Eq. (3.32) to 
obtain A. Indeed, Eq. (3.32) provides a good fit for a;„s where n = 2,3,4,5 and the 
resulting ojq is close to the expected asymptotic value 2(1 — i)7rT. To analyze separate 
contributions to the conductivity, we split the R.H.S of the sum rule (1.4) into two terms: 

si = -KIm{uJo), and S2 = -2KImA + 2 ^ Im (r^ - Kooq) . (4.1) 

n=l 

where in practice we set nmax = 5. In Fig. 2, we plot a, s, si, S2 and the total sum s = si+S2 
(normalized by 2ttK) versus T. We extract the conductivity via the analytic results of 
Ref. [35] (see also Ref. [33]): 

a = 27rA'Te^(^«). (4.2) 

Obviously from the plot, the R.H.S of the sum rule (1.4) s (dashed line), calculated 
numerically, is close to the conductivity cr(solid green line) given by Eq. (4.2) in the tem- 
perature range we are considering. We take that result as a numerical evidence that the 
sum rule (1.4) applies to the "soft-wall" model. Moreover, si is linear in T and has no c 
dependence. Thus, we could interpret si as the contribution to the conductivity from the 

^"With this choice, J/ip decay constant and the ip' mass and decay constant are ofT by nearly 20%. This 
is because while there is only one parameter c in the "soft-wall" model, the spectrum of the quarkonium 
at zero temperature is controlled by both heavy quark masses and the string tension (or Aqcd)- A more 
realistic holographic model of charmonium addressing this issue can be found in Ref. [33]. 
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thermal AdS background with no account of confinement effect introduced by parameter 
c. We also observe that S2 is always negative. That can be thought of as a reflection of the 
physical fact that the presence of bound states reduces the number of the charge carriers 
in medium and lowers the conductivity. That effect is quantified by ■52- Finally, we note 
from Fig. 2 that S2 in Eq. (4.1) is dominated by sa = —2Klm./\ term in the range of the 
temperature we are studying. That would mean that in some cases, one may be able to 
use —K\m. {loq + 2A) reasonable estimate of a. 

5 Summary and discussion 

We have shown that the current-current correlator in a theory with holographic dual de- 
scription can be represented as a convergent infinite sum over the quasinormal mode poles 
Eq. (1.1). We have established the convergence by deriving the asymtpotic behavior of the 
quasinormal mode frequencies and residues Eq. (3.29) using the WKB, or phase integral, 
approach. 

We have also established a sum rule relating conductivity a to the convergent infinite 
sum over quasinormal modes Eq. (1.4). We have checked this sum rule in the exactly 
solvable case of the A/" = 4 SUSY Yang- Mills theory. We studied the non-trivial example 
of the soft-wall holographic model numerically and found that the sum rule is in good 
agreement with analytically known value of the conductivity, and that the sum over the 
quasinormal modes is quickly saturated by a few lowest terms. 

5.1 Spectral function 

Using representation Eq. (1.1) for Gn we can also obtain a corresponding convergent rep- 
resentation for the spectral function: 



We have expressed the "gap" ojq and "offset" A parameters of the quasinormal modes 
in terms of the singular point of the Schrodinger quation ^oo and the corresponding Stokes 
constant S, Eq. (3.29), (3.30). Further insight into the significance of (or wq) and S (or 
A) may be obtained if one assumes that asymptotic expression of T){oj) in Eq. (3.28) can be 
continued to the real axis Arg(u;) = 0. Let us further assume that U{z) and h{z) are even 
functions of z. Then ^{z) defined by Eq. (3.7) is an odd function of z. Consequently, U{^) 
is an even function of ^. One can then argue that the corrections to Eq. (3.13) should be in 
even powers of uj~^ . However, as p{u}) is an odd function of w, those power corrections may 
not affect the asymptotic behavior of p(a;)^^. We could then use Eq. (3.13) and Eq. (3.28) 
to study the asymptotics of the spectral density: 



^^This is in agreement with the resuhs of Ref. [36] that spectral densities have no power law corrections 
in asymptotic expansion if the OPE of Euclidean correlators are free from non-analytic terms. 




(5.1) 



p{uj) ttKuj^ 1 + 2Im {Se^"^^°° ) 



ttKu^ 1 + 2e~^^^'lm (^e^^'^^^) 



(5.2) 
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where ^oo = ^B. + The first term in the square brackets, 1, on R.H.S of Eq. (5.2) 
is expected as p{uj) will asymptotically approach zero temperature limit. The next term 
explains the observation made on the basis of the numerical studies of Ref. [19] that "fi- 
nite temperature result oscillates around the zero temperature result with exponentially 
decreasing amplitude." The author of Ref. [19] argues that such behavior is intimately con- 
nected with the analytic structure stemming from the quasi-normal modes. Indeed, since 
Imwo < and thus > 0, Eq. (5.2) shows that 2^ij and 2^/ correspond to the oscillation 
frequency and the damping rate respectively. Our analysis suggests that such phenomenon 
is quite generic for theories with a gravity dual. 

One can easily check the correctness of Eq. (5.2) with Eq. (2.12) for the A/" = 4 SYM 
in the strong coupling limit where the Green's function is known analytically [24]. In 
addition, one can extend our analysis to other channels, e.g., the shear channel, as well. 
For example, for AT = 4 SYM in the strong coupling limit, again, ^oo = (1 + ^)/4T, we 
then predict the damping rate of the corresponding spectral density to be 1/(2T) while by 
fitting numerics, the authors of Ref. [37] obtained a damping rate of .46/T. 

In passing, we also note that due to the asymptotic behavior in Eq. (5.2), the integral 
over uj'^^~^ 6 p{u:) , where 6p{uj) = p{uj,T) — p{uj,T = 0), is convergent for any positive 
integer n. This suggests that for theories with a gravity dual, one could establish a family 
of f-sum rules [38] as discussed in Appendix C. 

5.2 Conductivity 

An insight into the meaning of the conductivity sum rule can be obtained by assuming a 
naive representation of the Green's function G/j(a;) in terms of the quasinormal modes: 



This representation ignores the fact the the sum is divergent. In a certain sense, the 
convergent representation (1.1) is a regularized version of the naive representation (which 
is indicated by letter R in (5.3)). Taking imaginary part in Eq. (5.3) and using Kubo 
formula (2.3) we would find 



Again, this sum rule ignores divergence of the sum. We can think of Eq. (1.4) as the 
regularized form of this naive sum rule. 

One could interpret the naive sum rule (5.4) as an expression of the following physical 
picture. Consider the behavior of quarkonia-like resonances as a function of temperature. 
As the temperature is increased the resonance poles in the Green's function move into (the 
lower half of) the complex plane. The residues r„, starting off as the real decay constants 
at T = 0, acquire their imaginary parts at finite temperature and are thus related to the 
process of "melting" or dissociation of the resonances. As a bound state dissociates, the 
conductivity receives contribution from the freed charge carriers. 

Although this picture is intuitive, its usefulness is limited, as the actual example we 
considered in Section 4 shows. We find that the dominant contribution to the conductivity 




(5.3) 




(5.4) 



n 
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comes from the first term — i^Im(a;o + 2A) of the sum rule Eq. (1.4). This term could be 
thought of as the combined contribution of the whole tower of resonances as it is a result 
of the regularization of the divergent sum in Eq. (5.4). 
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A Redundant poles and u = —2in7rT 

We now discuss a subtle issue in the context of gauge/gravity dual on how to define Gr{oj) 
at following points: 

w = -2m7rr, n = l,2, .... (A.l) 
We consider the Frobenius power series expansion Eq. (A. 2) near zh'- 

oo 

V±{zM = {z- zhY^ - ^^)' (^-2) 

j=0 

where the indicial exponents r± = zizioj/iirT and [39]: 



Here, Fj is a linear combination of Cj_^, ■ ■ ■ ,c^ 





Fr = Y^[{k + r_)aj_fc + uj^^j-k]Ck (A.4) 

k=0 

and aj , 13 j have no uj dependence [39] : 

,/ oo , n2 o° 

{Z - Zh){B' +-^) = Y1 «'^^(^ - ^^)'' = E - ^Hf- (^-5) 

A;=0 fc=0 

As a result of Eq. (A. 3, A. 4, A. 5), the coefficient c~{uj) will be a meromorphic function of 
UJ with simple poles at uj given by Eq. (A.l) for j > n. However, those singularities can 
be cured naturally by suitably choosing the overall constant Cq . For example, one may 
define^^: ^ 

= r(i - iu/{2TiT)) ' ^^-^^ 



This trick has been used to analytically continue the Gauss hypergeometric function in its parameter 
space. 
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then the Probenius solutions (A. 2) are regular in the entire complex u plane. Consequently, 
those points listed by Eq. (A.l) will not, in general, lead to any additional singularities of 
Gii{uj). Noting from Eq. (3.6) that Gji{u}) has no dependence on the overall normalization 
of the solution V-{u},z), a different choice of Cq will not affect resulting Gfi{uj). 

In fact, those special points have been known as "redundant zeros (poles)" [40] since 
long time ago in the context of the non-relativistic scattering. It has been shown [41, 42] 
that for the Schrodinger potential with the exponential tail: 

U{C) ~ e-y^ C^oo (A.7) 

the infalling wave solutions ^{^,uj) will have simple poles at u = —my/2. One can check 
from Eq. (3.10) and Eq. (3.7) that in our case, y = iirT, thus again we have Eq. (A.l). 
Historically, those poles are called "redundant poles" because they do not represent the true 
resonant states. In the context of the holographic correspondence, the name "redundant 
poles" may still be appropriate as they are not related to the singular points of the retarded 
Greens function Gji{u}). As a result, the singularities of Gfi{uj) are only due to simple poles 

at iOnioJn)- 

B The Stokes constant for a regular singular point 

In this section, we will derive the Stokes constant with respect to a regular singular point 
of a Schrodinger-type equation 

+ (a2 _ Uo{y)) ^{y) = 0. (B.l) 
Without losing generality, we set y = to be that regular singular point such that: 

limy^Uo{y) = l^ (B.2) 

and A to be a real positive large parameter. Eq. (B.l) has two Frobenius series solutions: 

i^iiy) = y^^ gi{y) (B.3a) 
iJ2{y) = y^' g2{y) (B.3b) 

where 5i(y), 52(2/) are power series of y. The indicial exponents are the roots of the equation: 

f-f-f = 0. (B.4) 

One may note from Weda's theorem that 

/+ + /- = 1. (B.5) 

A solution of Eq. (B.l), ^(y), can be expressed as a linear combination of ^i{y),^2{y), 
i.e., 

^(y) = ci*i(y) + C2^'2(y). (B.6) 
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Figure 3. The schematic plot of the Stokes diagram for Eq. (B.l). One can compare it with 
Fig. 1(b). In large |a;|, two turning points ^t,^t shrunk to ^txj- What is more, the distance 
between ASi and ^5*2 (or AS2 and AS[) also vanishes. Therefore, schematically, we only draw two 
anti-Stokes lines in blue and two Stokes lines in red dashed line in this figure. We also plot the path 
r in black thick line and F' black thick dashed line. 

From Eq. (B.3), we also have: 

^(ye^^^) = cie^^^-^+'f i(y) + cse^^^^- ^2(y); (B.7a) 

= cie-2^^^+^i(y) + C2e-2^^-^-^2(y). (B.7b) 

Multiplying Eq. (B.7a) and Eq. (B.7b) by e-''*(-/'+ and e''*(^++^-) respectively, then 
adding the results together, we have a connecting relation [43]: 

2cos[7r(/+-/_)]^(y) = e-"^(-^++^-)$ (ye2"*)+e"*(^++^-)^(2/e-2-*) = _[$(ye2-^)+^(ye"2-*)]. 

(B.8) 

If we solve the Schrodinger equation (B.l) using WKB approximation, we find two turning 
points at y = zizl/X. For A — )• 00 these points approach the singularity at y = 0. In 
the region of y not very close to the origin, where ^ \U{y)\, the solution ^'(y) can be 
approximated by a linear combination of two WKB solutions: e^*'^^. 

We plot the Stokes diagram schematically in Fig. 3. For large A the region where the 
WKB approximation breaks down shrinks to the origin. This region includes both turning 
points and the singularity at y = 0. Since we are working outside that region, so that 
A^ ^ \U{y)\, we can represent this non-WKB region by a single point at y = 0. Only four 
anti-Stokes lines emanate from this region as shown in Fig. 3: two from each turning point, 
following the real axis in positive and negative directions. 
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To determine the Stokes constant, we will consider a WKB solution defined by its 
value along the real axis: 

qfQ(y) = e'^y, whenArgy = 0. (B.9) 

When continued counterclockwise from the positive real axis to the negative axis along 
r', ^'o is unchanged as e*^^ is subdominant compared to e~*^^ in the upper half plane. If 
we go on continuing ^o{y) along T from the negative real axis to the positive real axis, we 
will have: 

^oiye'^""') = e'^y + Se-'^y. (B.IO) 

due to the Stokes phenomenon. Similar, when ^o{y) is continued clockwise from the 
positive real axis to the negative real axis along F, we have: 

$o(ye-^^) = e'^y - Se-'^y. (B.ll) 

Furthermore, when ^'o(ye~'^*) is continued clockwise from the negative real axis to the 
positive real axis along F', we obtain: 

^o(ye-2"^) = (1 + S^)e'^y - Se-'^y. (B.12) 

Substituting Eq. (B.IO) and Eq. (B.12) in Eq. (B.8) and comparing the coefficient of e^*^^, 
we have: 



±2z cos 



(B.13) 



the desired Stokes constant. 

As pointed out in Ref. [43], if y is an irregular singular point, in general, it is still true 
that there exist solutions with the properties: 

V^i(ye2-) = e2-^^+^i(y) 

Mye^n = e^'''^-My) (B.Ub) 

with f± called "circuit exponents" [43] of the singularity, in analogy with indicial exponents. 
One then observes immediately the connecting relation (B.8) is also true for y being an 
irregular singular point. Because of that, one may generalize the present approach to 
determine S to the cases that y is irregular as well. 

C f-sum rules from holography 

We now derive a family of f-sum rules [38] for theories with a gravity dual. Due to the 
representation (1.1), we have the following dispersion relation: 



6Gni^u.E)=nu:)+ H (c.l) 

y_oo 27r W - lUJE 



where V{uo) is a polynomial of uj. We denote by "J" the difference between the value 
of a function at temperature T and the value of that function at zero temperature. For 
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example, 6Gfi{uj,T) = Gji{uj,T) — Gji{uj,T = 0). Applying the Borel transformation to 
Eq. (C.l), we have: 

B^,t^5GR{ojE) = -2tB —5p{uj)s\n{ojtB) (C.2) 

where we have used the fact that p{ijj) is an odd function of u. We then consider the 
asymptotic expansion of 5Gii{iojE) 

5GR{iujE) = when oo (C.3) 

n=0 ^ 

where the coefficients 5hn may be calculated from OPE. Applying the Borel transformation 
to Eq. (C.3), we have: 

n=0 ^ 

Therefore 5hn can be extracted by comparing the Taylor expansion coefficients of Eq. (C.4) 
with that of Eq. (C.2): 



6hn = lim 



—dp[uj)sm[u;tB) 

0+ 



'■oo J 

2(-l)« / |-a.2-i<5p(a;). (C.5) 



As the integral is convergent, we have interchanged the sequence of taking the limit and 
the integration. In literature, Eq. (C.5) is related to the f-sum rule. From the definition of 
p{t) = {[J^{t), J*(0)]) and Heinseberg's equation of motion, we have: 

^^p{t) = (-^)^"'^([- • • [J'jt), J-(0)], T°°(t)], . ,r°°(t)] ). (C.6) 

2n 2n-l 

By the Fourier transformation and taking f — )• limit, we have: 

oo J 

^u'--^p{u) = {[...[ J\0), J^O)], roo(o)], . . . , too(o)]). (C.7) 

2n 2n~l 

Comparing sum rules (C.7, C.5), we have established an expression of 6hn- 

6K+1 = {-IT5{[. . . [ J'(0), J^(0)], T°°(0)], . . . , r°°(0)] )r. (C.8) 

2n 2n-l 

Finally, we list all the formulas of the Borel transformation that we have used 

Byt, (^) = tse-^'-, B,/,^ (3r) = ^lA. (log-) = -1 (C.9) 

for reference. We have also used the property that t3 gives zero when acting on polynomials 

of bJ. 
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